Fermions in Brans-Dicke cosmology 
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Using the Brans-Dicke theory of gravitation we put under investigation a hypothetical universe 
filled with a fermionic field (with a self interaction potential) and a matter constituent ruled by a 
barotropic equation of state. It is shown that the fermionic field (in combination with the Brans- 
Dicke scalar field <p{t)) could be responsible for a final accelerated era, after an initial matter 
dominated period. 
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I. INTRODUCTION 

Einstein's General Relativity (GR) is the usually in- 
voked theory for describing the evolution of the funda- 
mental space-time variables in cosmological models p]]. 
As an alternative to GR we have the Brans-Dicke the- 
ory of gravitation [l], a scalar-tensor formulation that 
rules the gravitational phenomena through the interplay 
between the metric tensor and a scalar field tp that con- 
trols the intensity of the gravitational constant G. On the 
other hand, cosmological models that include dark energy 
are strong candidates for explaining thepresent regime 
of positive acceleration of the universe [l], Going in 
that direction one possibility is to consider fermionic 
fields as gravitational sources for those accelerated uni- 
verses; these sources has been investigated using sev- 
eral approaches, with results including exact solutions, 
anisotropy-to-isotropy scenarios and cyclic cosmologies 
(see, for example [j, Hj-0])- Recently, these authors Q 
proposed a cosmological model, based on GR, in a dissi- 
pative Universe and showed that in a young universe sce- 
nario the fermionic field produces a fast expansion where 
matter (included via a barotropic equation of state) is 
created till it starts to predominate and the initial ac- 
celerated period gives place to a decelerated era. In this 
case the fermionic field plays the role of the infiaton in 
an early period and the role of dark energy for an old 
universe; without the need of a cosmological constant 
term or a scalar field Q- In an exclusive old universe 
scenario an initially matter dominated period gradually 
turns into a dark (fermionic) energy period when an ac- 
celerated regime rules for the rest of the system evolu- 
tion. One important point is that it was considered a 
self-interaction term for the fermionic constituent, in the 
form of a function that included the Nambu- Jona-Lasinio 
potential as a particular case [8]). Testing those proper- 
ties in a universe ruled by Brans-Dicke gravitation would 
be of interest and this is the main focus of this work. 
The manuscript is structured as follows. In section II we 
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make a brief review of the tetrad formalism used to in- 
clude fermionic fields in a dynamical curved space-time. 
The Brans-Dicke field equations for an isotropic, homo- 
geneous and spatially flat universe are derived in section 
III where we present the analysis of the different sce- 
narios in which the fermionic constituent would, in prin- 
ciple, answer for the transition to accelerated periods. 
Finally we display our conclusions. The metric signature 
used is (+,—,—,—) and units have been chosen so that 
8ttG = c = ^= 1. 



II. 



FERMIONS AND GRAVITATION 



When one tries to include fermions in a gravita- 
tional model what must be taken into account is that 
the gauge group of general relativity does not admit a 
spinor representation. The tetrad formalism solves the 
problemQ, 0, ©, Ejj. Following the general covariance 
principle, the metric tensor satisfies 



3^7 



0,1,2,3 



(1) 



where denotes the tetrad or vierbein and r\ a h is the 
Minkowski metric tensor; latin indices refer to a local 
inertial frame whereas Greek indices to the general co- 
ordinates system. The next step to couple the fermionic 
field to gravity is to construct an action for the model; 
we start reminding that the Dirac lagrangian density in 
Minkowski space-time is 



(2) 



where m is the fermionic mass, V = V^7° denotes the 
adjoint spinor field. V is a function of ip and ip, and 
represents a fermionic self-interaction. The general co- 
variance principle imposes that the Dirac-Pauli matrices 
7° must be replaced by their generalized counterparts 
r M = e^7 a , where the new matrices satisfy the extended 
Clifford algebra, {r^,r i/ } = 2g» v . In a second step we 
need to substitute the ordinary derivatives by their co- 
variant versions 

(3) 
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where the spin connection is given by 
1 



<T A 



(4) 



with rj^ A denoting the ChristofFel symbol; the generally 
covariant Dirac lagrangian then becomes 

L D = ^®T»D^- (D^)T^] - m^/- - V. (5) 

The field equations are then obtained from the total 
action 



l—g Lj< d x, 



(6) 



where Lt = Lq+Ld+Lm is the total lagrangian density. 
Lm is the lagrangian density of the matter field and Lq 
is the free gravitational lagrangian density. In the Brans- 
Dicke theory [l| we have 



{Re av -a 2 oje av (W(p) 2 }, 



(7) 



where R is the curvature scalar and ip is an scalar 
field that controls the gravitational constant (Gn) 
intensity [lj]; a and ui are constants[l|, Ld is the Dirac 
lagrangian density ([5]) . Using the variational principle we 
obtain then the Dirac equations for the spinor field (and 
its adjoint) coupled to the gravitational field 

iT^D^-m^^^L = 0, zZV0P"+m^+^ = 0. (8) 
dip dip 

Analogously, the Brans-Dicke gravitational dynamics 
emerges from the total action ([6]) 



G ^ v 



(9) 



where 



G^ v = e aH> {R iiV - -g^R ~ ^(g^Dip - V„V^) 



(10) 



and is the energy-momentum tensor of the sources of 
the gravitational field: T^ v = Tg'+Tjg. The symmetric 
form of the fermionic energy-momentum tensor is given 
by 

- D^Vtp - D^Vip] - g^L D . (11) 

These equations rule the dynamics of a Brans-Dicke uni- 
verse filled with a fermionic and matter sources. 



III. BRANS-DICKE AND DIRAC FIELD 
DYNAMICS 

Friedmann-Robertson- Walker (FRW) metric mirrors 
the homogeneity and isotropy properties of the universe. 
The space-time interval is usually written as 



ds 2 = dt 2 - a(tf(dx 2 + dy 2 + dz 2 ), 



(12) 



where a(t) is the cosmic scale factor. In the FRW metric 
the tetrad components (1) become 



y - — !— 

°' 4 " a(t) 



<5 M 



(13) 



and the Dirac matrices turn out to be 



r° = 7 °, 



r = 



a(t) 



r 5 = 



^ r r l r 2p3 = 

. (14) 

from which the spin connection (see eq. (4)) is ob- 
tained; the non-zero components are Sli = ^a(^)7 l 7°- For 
this isotropic and homogeneous universe the fermionic 
field becomes an exclusive function of time; therefore the 
Dirac equations (|HJ| for a non-massive fermionic field be- 
come 



i/j+-HiP + tl °-= 
2 dip 



0, V + ^~^7° = 0. (15) 



Besides the Dirac field, we have considered as grav- 
itational source a matter constituent that is ruled by 
a barotropic equation of state, namely p m — b m Pm, 
where b m is the barotropic coefficient (0 < b m < 1); 
then the total energy density p is given as a sum of 
these individual contributions, i.e., p = po + pm — 
Too, where T^ v is the total energy-momentum tensor of 
sources. The energy density p satisfies the conservation 
law p + 3d/a(p+p) = 0. Besides, the non-vanishing com- 
ponents of the fermionic energy-momentum tensor follow 
from (12) yielding 



(Td)° 



V 



PD, 



(To) l = V- — - - --= = -p D , 
dip 2 2 dip 



(16) 



(17) 



which are exclusive functions of ip and ip. Combining 
the Dirac and Einstein equations one can obtain an in- 
dependent conservation law for the energy density of the 
fermionic field Q. This implies into a decoupled conser- 
vation equation for the energy density of the matter con- 
stituent, i.e., 



Pn 



Pn 



0. 



(18) 



where H — a/a is the Hubble parameter. 

In order to analyze the cosmological solutions of our 
model we have to define first some sources properties. 
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For the fermionic field we consider a self-interaction po- 
tential V; this can be modeled as an exclusive function 
of the scalar invariant (ipip) 2 ; following the Pauli-Fierz 
theorem 0,(1] we have 



V = [(VV) 2 



(19) 



where n is a constant real number. A particular case of 
the Nambu-Jona-Lasinio potential H is obtained when 
n = 1. In fact, the fermionic field behavior can be clas- 
sified according to the value of the exponent n. In par- 
ticular, for n < 1/2 the pressure of the fermion field 
is negative and it could represent (in a universe ruled 
by Einstein gravity) either the infiaton or a dark energy 
constituent I7|. The fermionic energy density is given by 
Pd = [(V'V') 2 ] • O n the other hand the second con- 
stituent, the matter field, is described, as mentioned in 
the precedent section, following a barotropic equation of 
state[7|. After some algebraic manipulation we can put 
the model dynamics in the following form: 



2- + H 2 + c?(p + 2aHip H ip 2 

a 2 



a (p + 'daHp 



3 + 2w 
3(2n- 1)(^V) 2 ™ 



V + p, 
- 3p, 



(20) 



p rn + 3H(p m +p m ) = 0, 



<Z dtp 



ip + -Hip 



(21) 



(22) 



z^7° = «23) 



dip 



We consider a pressureless matter field (p m — 0), so 
that one can obtain from equation (|22p that p rn (t) = 
Pm(0) / a(t) 3 ■ Furthermore, the equations (f2"3"f can also 
be integrated for the potential (fT9f resulting that the 
fermionic bilinear evolution is governed by [iprp](t) = 
[tptp](0)/a{t) 3 . The remaining equations {11 and ([2Tjl 
constitute a highly non-linear system of differential equa- 
tions and we proceed to solve it numerically. We analyze 
the time evolution of our model choosing first the condi- 
tions for t = 0: 

a(0) = l, <x(0) = l, $0(0) = 0.001, p m (0) = l, 
93(0) = 1, <£(0) = 0.001. 

These conditions characterize qualitatively an initial pro- 
portion between the constituents; an era when matter 
predominates over the fermionic density. Besides that, 
we have to specify the magnitude of the remaining pa- 
rameters: we suppose initially that a — 1.0, w = 4x 10 4 
[ll| and a value for the potential power n, n = 0.2 Q. 

These choices are reference values that permit final ad- 
justments to follow several cosmological constraints, like 
the spectrum of the Brans-Dicke coupling w[ll[ and the 



present value of the Brans-Dicke scalar field ip. These pa- 
rameters can be in fact adjusted due to invariance prop- 
erties of the Brans-Dicke gravitational field equations, by 
using the following change of variables: a — > a, p — > 
(p, ap — ~f/tp + a where exp(— 7) = xq is associated 
to the asymptotic value of (po — (2lu + 4)/(2w + 3). 
These transformations show how a value of the Brans- 
Dicke parameter uj is linked to the definition of a new 
Brans-Dicke scalar field (p. In fact, after numerical in- 
tegration, it is possible to verify that these features 
are included in the ip evolution (see figure 1) that im- 
ply into the evolution of the gravitational "constant" as 
G(t) = (2w + 4)/[(2w + 3)£(t)]. 

In figure 2 it is plotted the acceleration field d as func- 
tion of time t for two different values of n. The results 
show that initially the universe is expanding with nega- 
tive acceleration, a period where the matter constituent 
predominates over the fermionic field. With the evolution 
of time we have increasing values in energy transference 
to the fermionic field (this is not happening directly but 
via the gravitational field a(t) and the scalar field <p(t), 
as the equations of motion (20-23) show). The negative 
pressure of ip(t) help in fact to promote a final acceler- 
ated period (also shown in figure 1) indicating that in 
this model the dark energy role would be played by a 
combination of the fermionic constituent with the scalar 
p(t). Other interesting results appear when we choose 
the power n to be in the neighborhood of values n w 0.33. 
In fact, for a fixed value of <^(0) and increasing values of 
n what emerges is a universe that fails to show a final 
accelerated period, even when the fermionic field still ex- 
hibits a negative pressure {po = (2n — 1) [(V'V') 2 ] )• On 
the other hand, for values n < 0.33 we will find a universe 
that is permanently in accelerated expansion. Another 
important remark here is that this general qualitative be- 
havior depends strongly on the initial value of the time 
derivative of the scalar field tp(t): what we verify is that 
increasing values of <p(0) promote an earlier entrance on 
the accelerated period. 

In figure 3 it is shown the behavior of the energy densi- 
ties of the fermionic field pd and matter p m as functions 
of time t. The numerical results show that eventually 
the energy density of the femionic field overcomes the 
energy density of the matter field, although this does not 
coincide with the instant when the universe goes into an 
accelerated period (opposed to what occurs in some Ein- 
stein gravity based models [H|)- Another feature showed 
by the numerical results is that for larger values of n it fol- 
lows that: (a) the energy density of fermionic field grows 
more slowly causing a larger decelerated period and (b) 
the energy density of the matter field has a less signifi- 
cant decay and also promotes a larger decelerated period. 
On the other hand, as figure 3 shows, both densities have 
decreasing values in time due to the permanent expan- 
sion of the universe (this was verified with plots of the 
scale factor against time). Again, all cases are strongly 
dependent on the exponent n of the self-interacting po- 
tential and one can obtain different behaviors, that are 
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in tune with the acceleration patterns presented above. 
Finally we verify that the scalar field ip approaches a final 
constant value that validates the accord between Brans- 
Dicke and Einstein gravitation for large t. 

As concluding remarks we stress that we have investi- 
gated the role of a fermionic field - with a self-interacting 
potential - in an old Brans-Dicke universe. We have 



shown that the fermionic field, in combination with the 
scalar (p behave like a dark energy constituent, promot- 
ing a final accelerated period after an initial era domi- 
nated by matter. These results depend strongly on the 
initial values of <p and on the power n present in the 
self-interactive fermionic potential. 
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